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I. INTRODUCTION. 



New data on the photoproduction and electroproduction of charm in deep inelastic scattering 
(DIS) recently obtained by collaborations of ZEUS and HI at the HERA accelerator, stimulate 
the discussion on the following problem: what is an interpretation of these data in the framework of 
perturbative QCD (pQCD)? At present the data show the experimental confirmation of qualitative 
predictions of pQCD. However, essential quantitative deviations of measured cross sections from 
the theoretical expectations were observed in some kinematical regions, i.e. the discrepancies in the 
form of Z}*-meson spectra were found. 

Exploring the perturbative theory in the processes under consideration is strictly justified due 
to quite large transverse momenta of D*-mesons, px > "^c- To the same moment, we cannot com- 
pletely avoid nonperturbative terms in the calculations, and the hadronization process of c-quarks 
produced in hard collisions of initial partons is usually described by means of a fragmentation func- 
tion (FF) D{z,fi), where the parameter fi determines the scale of factorization for the probabilities 
of perturbative subprocess and binding the quarks into the hadron. The hadronization of quarks is 
determined by the dynamics of confinement^. The information on the nonperturbative FF is taken 
from the data on the production of charmed particles in e"'"e~-annihilation The dependence of 
FF on /i is determined by the leading logarithmic approximation of QCD (LO) and corrections to it 
(NLO). Calculating the cross sections for the photoproduction of heavy quarks in the 0(a;Q;^)-order, 
some authors neglect the D{z,^) changes versus ~ p-r in the hadronization [Q. This dependence 
of FF was taken into account in other papers 0. In this way, the authors used the initial FF at 
/i = /io ~ 1 GeV in the forms of both the ansatze based on the reciprocity relation 0] and FF by 
Peterson et al. 0. 



*Another region, wherein phenomenological parameterizations are necessary because of the nonperturba- 
tive character of interactions, is the description of parton distributions in initial states, though the evolution 
versus the scale of factorization for the structure functions allows the consideration in the framework of 
pQCD. 
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In pQCD two approaches are used for the calculation of cc-pair production up to the 0{aa1)- 
terms. In the first, authors assume that the light quarks and gluons are the only partons present 
in the photon and proton, when the c-quarks are produced in the interactions of those partons. In 
this way the nonzero mass of c-quark is taken into account [^. 

The second approach suggests that the c-quark is the additional "active" flavor supposed as 
massless, and the c-quark distributions in the initial particles are introduced in the form of corre- 
sponding structure functions 

Note, that in the both ways the hadronization of c-quarks is described in the framework of 
fragmentation model, which results in the convolution of perturbative distribution of charmed quarks 
with the fragmentation function: 



dan* /"^ dcTccikx^lj) 



^^^^^:^dz, (1) 



dpT j2pt/v^ dkr 

where Dc^[)*{z, yu) is the FF normalized by the probability of c-quark transition into the D*-meson 
w{c D*), which is measured in the e+e~-annihilation [|| {w{c — > D*) = 0.22 ± 0.014 ± 0.014), 
and fi is the factorization scale for the perturbative partonic cross section d&cc/dkT- 

According to the model of heavy quark fragmentation, the basic condition of its applicability is 
Pt ^ '^c oi nice ttld providing a jet-like shape of measured events. However, the dominant term 
of statistics is integrated in the region of px ~ ttLc, where the model of (|I]) may not be explored. 

In the model offered below, we describe the production of cq quark state in pQCD and use the 
semilocal duality [TO] assuming that this quark state repeats the distributions of hadron with the 



accuracy up to the factor depending on the quantum numbers of cq and the hadron. In this way 
we consider the complete set of Feynman diagrams corresponding to the production of cg-system in 
the fixed order over a and as (so the leading order is the fourth one). Then we quite successfully 
describe the production of D*-mesons at pt ~ ttLc as well as reliably and correctly approach the 
results of fragmentation model at px ^ fric- Such the agreement of model estimates with the 
experimental data can be reached, if we add the term of octet-color state for the production of cq 
pair along with the contribution of the singlet-color state. 

In the framework of pQCD taking into account the full Q^- dependence of the matrix element. 



exact calculations give a good description of data on the DIS production of charm. This estimate 
is carried out up to NLO in refs. W^. Despite of the qualitative success in such the approach, 



quantitative deviations between the predictions and measured differential cross sections take place 
as in the photoproduction of charm. In present paper we calculate the differential cross sections for 
the D* production in DIS in the framework of model taking into account the mentioned contribution 
of octet operator into the production of charmed quark-light quark pair. 

II. FRAGMENTATION FUNCTION. 

Among the well-known parameterizations of FF, the model by Peterson et al. [§] argues the 
following form: 

V 2 {I — Z}' 

where is the normalization factor, and e is the free phenomenological parameter dependent of 
fi. This FF quite accurately describes the data on the production of B- and D-mesons in e~^e~- 
annihilation at high energies, where the fragmentation mechanism dominates JT2|. 



According to motivations by the authors of parameterization (H), the behaviour versus z in 
Eq. (H) is determined by the virtuality in the propagator of heavy c-quark. Indeed, let Pd* and Pjet 
denote the momenta of Z)*-meson and an associated jet, respectively. Then the denominator of 
perturbative propagator has the form 

Expanding this expression in small parameters and and defining z = ^ ^'l* — , we get 
- [Pd* + Pjet) ~ "^c ~ "^c ~ ^ 



1 — z ^ z 1 — z z ml 1 — z 

In the parameterization of (Q) the quantities N and e are treated as phenomenological and nonper- 
turbative parameters, though we can clear see that e corresponds to the mass ratio squared for the 
jet mass over the quark mass, whereas the jet mass logarithmically grows with the increase of total 
energy in the e"'"e~-annihilation. 



The consistent consideration of heavy quark fragmentation \T2\ with account for both the quan- 
tum numbers of meson in the final state and the QCD-structure of vertex on Fig. 1, leads to the 
following analytical expressiorQ: 

^--<^> = 274' (1- (1-4)' " ' ' ^ ' ^ ' 

2(1 - r)(4 - r + 2r^)z^ + (1 - 0^(3 - 2r + 2r^)z% (3) 

where r = niq/ {rriq + rric) and 



which corresponds to the square of wave function at the origin for two quarks in the framework 
of nonrelativistic potential model: (0(|))|ArR = 1^^(0)1^. The effective mass of light quark niq 
determines the form as well as the normalization of FF: 

Jo mi 
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where 



^The quark-meson vertex contains the form factor which makes the loop calculations free off any diver- 
gences: ultraviolet or infrared. So, the calculations are free off any regularization procedure, and we do 
not need any discussion on the regularization. This is very well known from ref. [l^. To the leading ap- 
proximation the quarks entering the meson have the same velocities and are on mass-shells (see ref. [T^). 
This takes place if we calculate the fast or valence components of the meson. The perturbation theory is 
useless, if one calculates the soft quark-gluon sea, which is out off current study, since we restrict the gluon 
virtualities and consider the fast degrees of freedom in the meson. According to the factorization theorem 
in perturbative QCD at high px the restricted number of diagrams contributes. It is well known that the 
only hard diagram (Fig. 1) contributes to the fragmentation in the special axial gauge (see ref. In 
the covariant gauge two hard diagrams do. The calculations of fragmentation functions were done in both 
gauges, and the FF does not depend on the gauge. 
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r(7 - 4r - Sr^ + lOr^ + 2r^) , , , 
+ ^T^e Hr) 



(5) 



In the model under consideration we put the probabihty of fragmentation w{c — > D*) independent of 
the scale /i/j, so that w remaines the constant value extracted from the experimental data mentioned. 

We argue the following motivation: The distributions of partons inside the heavy-light meson 
have distinct terms corresponding to the valence quarks. So, for the neutral charmed meson cu, 
say, in the infinite momentum frame these contributions are given by the functions depending on 
the fractions of longitudinal momentum and the transverse momenta of partons. The valence terms 
are determined by 

Plix,p±) =Pcix,p±) -Pcix,p±), (6) 
pUx,p±) =Pu{x,p±) -Pn{x,p±), (7) 

so that the functions yield the average fractions 

£p^dxx-p:{x,p^)^j^, (8) 

r A 
{xD = J (fp^dx x-pI{x,p±) ^ — , (9) 

with (xj!) + (x^) ~ 1, and A denotes the binding energy of heavy quark inside the meson. Then, in 
the model we neglect the dispersions of quark fractions and put the valence quark velocities equal 
to each other: Vc = Vu- Generally, we can write down the following representation for the Fock state 
of charmed meson: 

\D) = Wn.\icum + W^s)\icu)^s)), (10) 



where the quark pairs (cm) carry the mentioned fractions of meson momentum, and they belong to 
two possible color states: the singlet or octet, and the operators ^^(1,8) generate the quark-gluon 
sea inside the meson at the characteristic scale of bound state fiR. There are several qualitative 
features of this representation. First, the operators W are nonperturbative. Therefore the scale of 
interactions in W is about the scale of confinement Aqcd- In contrast to the heavy-light meson, 
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for the heavy quarkonium there are the counting rules due to the small relative velocity of heavy 
quarks, i.e. due to the smallness of typical inverse size of bound state with respect to the quark 
mass. In the heavy quarkonium the W(^s) contribution is suppressed by a power of relative velocity 



T3|| . In the heavy meson under consideration, the octet contribution is not suppressed, and it is 
of the same order as the singlet one. Second feature concerns for the excitations in the system of 
valence quarks. For the heavy quarkonium the colored P-wave states contribute, since their sizes 
are still small enough to propagate with no decay into the pairs of heavy mesons. In the heavy-light 
system the P-wave states are suppressed because they are heavy and have large sizes leading to an 
instability. Then, for the production of charmed mesons the contribution of singlet and octet S-wave 
state is essential only. The operators W^(i,8) are infrared and cannot be calculated in perturbative 
QCD. However, they do not influence onto the distributions of meson at the transverse momenta 
higher that the charmed quark mass, px ~ ttLc- To get rid of the nonperturbative quark-gluon sea 
from the calculations, we have to restrict the gluon virtualities kg > so that ^^qcdI lAi ^ 1- We 
suppose yU/j ~ m-p, and consider the massive light quarks with uiq ~ wip/2. 

Further, along with the production of color-singlet state of (cg)-system we have to take into 
account the production of color-octet state with the corresponding transition into the hadron {D*- 
meson). As we know from the experience on the consideration of production for the mesons contain- 
ing two heavy quarks, the octet term is suppressed by the fourth power of relative velocity v for the 
quarks inside the meson [|1^]. In the problem studied the velocity is not small and the contribution 
by the octet can be comparable with the singlet term. Due to the specifics of color amplitude for 
the e"'"e~ -annihilation the ratio between the production of color-singlet and color-octet (cg)-states 
is fixed, so that introducing the octet term for the fragmentation we have to substitute for {O^) 
in Eqs.(3) and (D by {O^ff) as follows: 

(O^^0 = (O(i)) + ^(O(s)), (11) 

where 

(0(8)) = ^ [-9^" + ^) {D\p)\{c^,X'^q){q^A'c)\D\p)) 
As we have already mentioned, for the mesons containing the light relativistic quark we expect 



(0(8))/(0(^) ~ 1, since f ~ 1. Then from Eq.(|TT]) we see that at such the ratio of operators the 
singlet production dominates in e^e^-annihilation. The form of differential cross sections is one and 
the same for the singlet and octet mechanisms in the fragmentation model. Therefore, we cannot 
distinguish the singlet and octet terms in the production of D*-mesons in the e"'"e~-annihilation. 
However, in deep inelastic scattering the situation is essentially changed in the region of pt ~ rUc, 
where the fragmentation is not the only mechanism, determining the production of heavy meson. 

Thus, we have quite a definite picture for the description of charmed meson production in 
perturbative QCD with the factorization of terms responsible for the contributions by the singlet 
and octet operators. 

Concerning the form of FF in this approach of factorization for the hard perturbative amplitude 
and nonperturbative binding operators, we note that the virtuality of heavy quark is equal to 

ml - {Po^ + P,.,f ^ ml. ((1 -rf-\- = (1 - (1 - r)z)\ 

so that ^ e. At an appropriate choice of mq and m^ we can reach the FF behaviour similar with 
that of FF by Peterson et al. in @). The logarithmic dependence of e mentioned above can be 
reasonably neglected in the region, where In— ~ 1, i.e. it is not large. 

III. PHOTOPRODUCTION AND ELECTROPRODUCTION OF CHARM. 

As we have shown above, in the e+e~-annihilation we obtain the reliable description of cross 
sections for the production of charmed particles, if we use the perturbative calculations with a 
nonzero mass of light quark. In this way the calculations are reduced to the evaluation of single 
diagram in a special axial gauge (see Fig. 1) |]12[. In other processes the number of diagrams is 
essentially greater than the mentioned one. The complete set of such diagrams in the 0(Q;a^)-order 
is shown in Fig. 2 for the subprocess of 7* (7 D* + c + q. Among this set, wherein the initial photon 
is on mass-shell for the photoproduction and virtual for the deep inelastic scattering, we can identify 
the diagrams corresponding to the hard production of c-quarks up to 0{aas) with the consequent 
fragmentation of heavy quark in the 0(Q;^)-order (for example, graphs 16 and 19). These terms 
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dominate at pt ^ ttLc, which reproduces the predictions of fragmentation model by Eq.([^) in the 
kinematical region specified. 

Let us study the process 7(7 ^ D* + X at a rather high energy s ^ 4M|,, to prove the vahdity 
of our statement on the factorization. The large subprocess energy is required to reach the limit of 
Pt ^ nT'D, where the predictions of the model can be compared with the fragmentation regime in 
the form of (|lD. 

The calculation results for the py-distribution of the D*-meson production cross section in the 
framework of our model in comparison with the fragmentation prediction (|lD are presented in Fig. 3. 
In Eq. (|l|) the Born approximation for the cc-pair production and the fragmentation function in the 
form of (3) have been used. We can see that the factorization for the hard cc-pair production and 
the probability of c-quark fragmentation into D* for the photoproduction cross section takes place 
at Pt > 20 GeV for the singlet state and at pt > 40 GeV for the octet term at the photon-gluon 
energy equal to 200 GeV. We find that at the lower energy of photon-gluon interaction, say, at 
a/S = 40 GeV, the fragmentation regime works at pt > 12 GeV for the singlet state production, but 
the differential cross section for the production of octet state cannot be described by the factorized 
form of fragmentation since at such the energy the higher twists over pt dominate completely in 
the whole kinematical region. 

Thus, the model approaches the fragmentation regime at rather large transverse momenta. At 
Pt ~ rric the dominant contribution comes due to the diagrams of nonfragmentational kind. In 
these diagrams the independently produced cc and qq pairs after rescattering are transformed into 
hadrons. We call them as the recombination diagrams (see, for example, graphs 5, 6, 7 in Fig. 2). 
The fragmentation and recombination diagrams have different forms of asymptotic behavior at large 
Pt- So, the recombination contribution at Pt ^ ttLc drops like l/pj' in contrast to the fragmentation 
term decreasing as 1/p^- 

The fragmentation approach is valid at pt ^ ttiq and cannot be applied at pt ~ ttiq (see Fig. 3). 
To describe the region of pt ~ rnq we have to take into account subleading terms over 1/pt, i-e. 
the higher twists. The offered model allows us to get correct factorization at high pt and to account 
for the subleading terms on the basis of gauge invariance dictating the consideration of complete 
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set of diagrams to the given, fixed, order in ag- 
in the framework of fragmentation approach, the LO calculations must take into account the 
resolved photon contribution. In the model under consideration, which involves the higher orders, 
this contribution is also included by means of diagrams, where the splitting of photon into the 
pairs of light quarks is taken into account perturbatively at rather high pt, and, hence, at large 
virtualities|. 

An analysis of shown that the logarithmic corrections to the c-quark production become 
essential at pt > 20 GeV. Thus, we can believe that at moderate px'- In^ < 1, the fixed order 
calculations (LO, NLO) are valid. 

Remember, Eq.(0) cannot be applied at pr ~ rric because of two following reasons: 

1) The limit of pt ^ is not justified at px ~ fric, 

2) The higher twists over pt must be taken into account beyond the fragmentation regime leading 
to Eq.(g). 

In contrast, the offered model under consideration can be used in the kinematic region pt ~ rric 
too, where Eq.(0) does not work, and it reproduces (|l]) at pr ^ riic. 

Thus, exploring the gauge invariance to the fixed order, we extend the fragmentation term up 
to the complete set of diagrams to cover the higher twists over p^ ~ ttLc. 

We would like to stress that the perturbative calculations for the production of cg-system can 
reasonably describe the region of pr ~ ttLc, if we suppose the appropriate effective mass of light 
quark rriq. In the calculations we chose the value rriq = 0.3 GeV ~ mp/2. Then the gluon virtualities 
are greater than |A;^| ~ + m^, which allows us to use pQCD and to stay beyond the region of 
infrared divergency. 

Remember that at px > 12 GeV we normalize the absolute value of cross section due to relation 
(^, assuming that the probability of fragmentation does not depend on the scale. 

In contrast to e+e~-annihilation, in both the photoproduction and the electroproduction we find 



^We check that the nonperturbative splitting of photon into the pair of charmed quarks due to the 
contribution of vector bound states can be neglected in the kinematical region under consideration. 
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another ratio of fragment at ional and recombinational contributions of octet state in comparison 
with the ratio for the singlet. So, if (C'(8))/(0(|)) ~ 1, then the recombination contributes to the 
production of octet with the same order of magnitude as for the singlet production. Therefore, at 
the transverse momenta pt ~ rric, where the recombination dominates, the values of cross sections 
for the singlet and octet are close to each other, when at pt ^ rric, where the basic role is played 
by the fragmentation, the octet state term is suppressed as 1/8, according to Eg . (|TTD . 

We emphasize that the introduction of octet term does not practically change the cross section 
for the production of D*-mesons at px ^ itLc, however, it essentially contributes in the region of 
Pt - rric. 



A. The photoproduction of L'*-mesons. 

We fix the value of fragmentation probability w{c ^ D*) taken from the data on the D*-meson 
production in e"'"e~-annihilation at LEP, and we put the choice of fiR, rrig, rric, so that Eq.(^) 
definitely determines the value of {O^^^ {hr}) operator matrix element. So, if 

^J'R = rriD*, 
m„ = 0.3 GeV, 

(12) 

rric = 1-5 GeV, 
w{c^ D*) = 0.22. 

then 

{0'"{mD*)) = 0.25 GeVl 

Such the choice has allowed us to describe the data on the photoproduction of D*-mesons. 
The results of calculations are shown in Fig. 4 in comparison with the ZEUS data [|I|] measured 
in the following kinematical region: px > 2 GeV, -1.5 < r] < 1.5, 130 GeV < W < 280 GeV, 

< 1 GeV^, where W is the total energy of 7p-interactions, is the photon virtuality, and rj is 
the pseudorapidity of D*-meson. 

We see in Fig. 4 that the value of production cross section for the D*-mesons cannot be explained 
by the singlet term only. The relative value of the octet contribution (0(8))/ {O^) = 1.3 allows us 

11 



to reach the good agreement with the experimental spectra of D*-mesons. We put the attention to 
the mentioned distinct difference in the behaviour of distributions over the transverse momentum 
for the singlet and octet terms. The contribution of octet state is enforced in the forward direction 
[i] > 0). At the chosen values of operator matrix elements the contributions of singlet and octet are 
close to each other at pt ~ tUc, while at pt ^ rric the octet term becomes less essential than the 
singlet one. 

We can see that the accounting of the octet contribution essentially improves the form of pre- 
dicted cross section distribution over the pseudorapidity in the forward direction. Such the im- 
provement is due to the fact that the recombination diagrams dominant at pt ~ ttLc are enforced 
in the octet part of the cross section by the color factor. 

The interaction of virtual c-quark with the light valence quark in the initial hadron can be taken 
into account easily, as we will show below, but this process becomes essential only at low energies 
of charm production near the threshold. 

The scale of factorization for the gluon structure function in the proton has been fixed fip = 
2m D*- Calculating the differential cross sections we suppose two values for the scale of factorizing 
the matrix elements for the quark operators (0(i^8)): fiR = mu* (the upper line) and fiR = 2m£,* 
(the down curve). We see that the experimental data on the photoproduction of D*-mesons at 
HERA can be described in the model in the better way if we put fip = 2m d* and = m^*. 

The predictions for the photoproduction of D*-mesons in the low energy region studied by ZEUS 



T§: Pt>2 GeV, -1. < r/ < 1.5, 80 GeV < < 120 GeV, < 0.01 GeV^ are shown in Fig. 5 



In the calculations we use the same values of and (0(8)) as they stand above. 

Unfortunately, the present experimental data do not yet allow us to extract the value of (0(8)) 
with a low uncertainty, and with the mentioned variety in the factorization scales we can write down 
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(0(8)) ^ 0.33 - 0.49 GeV^ 

It is quite clear that changing the parameter (0(|))/(0(8)) we can enhance or decrease the con- 
tribution of higher twists in the region of px ~ In the calculations we neglect the difference 
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between the spectra of vector and pseudoscalar octet states, whose contributions are effectively 
summed up in the corresponding operator. Furthermore, we do not consider a possible influence 
of P-wave color-octet states, since their contribution is suppressed in comparison with the S-wave 



contributions, as we saw in similar studies ||16 



B. The photoproduction of Ds- and D*-mesons. 

In the model under consideration we predict the cross section for the photoproduction of Dg- 
and D*-meson. The general features of these processes were discussed in our work |jT^. A more 
detailed analysis is given below. 

The probability of c-quark fragmentation into Ds or D*, Ws{c Ds,D*) = 0.1, was fixed from 
the data on the D^-meson production in the e"'"e~-annihilation at LEP ||18|| . 

Using this value of Wg and replacing the effective valence quark mass rriq = 0.3 GeV by = 
0.5 GeV in the amplitude of D*-production, we calculate the expectation of differential cross section 
for D*. The analogous technique was used for the calculation of Dg production cross section. In 
the both cases the value 1.3 was chosen for the relative octet contribution such as in the D*- 
meson production. It is worth to mention that the flavor SU (3)-symmetry for the u, d, s quarks 
is broken. So, to obtain the experimentally measured fraction of charmed-stranged events in the 
charm production, Wg, we put the effective operator {O^^^ {mj:,*)) for the c-quark fragmentation into 
Dg to be approximately twice greater than {O"^^ ^ {m £,*)). 

In this work we present the differential cross section for the sum of Dg- and D*-mesons because 
the experimental discrimination between the events with Dg and D* is problematic. The reason is 
the necessary registration of a soft photon coming from the dominant mode of Dg D*'~f decay, 
which is quite difficult task. 

The results of our calculations in comparison with the ZEUS experimental data ||19| are presented 



in Fig. 6 for the following kinematical region: 130 < < 280 GeV, < 1 GeV^ pr > 3 GeV and 
-1.5 <T] < 1.5. 

We see that the predictions for the Dg and D* photoproduction cross section are in a good 
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agreement with the experiment. 



C. Electroproduction of Z?*-mesons. 

In contrast to the photopro duct ion, we have to take into account the exact dependence 
in DIS for the calculation of D* production. In this respect we evaluate the tensor of leptonic 
currents squared and averaged over the polarizations to contract it with the photon-gluon vertices 
in the charm production calculated numerically at arbitrary of photon. Again we consider the 
production of cg-pair in the 0(ag^ag)-order. 

The used technique of calculations allows us to evaluate the cross sections for the production 
of Z)*-mesons in the case when the initial photon is essentially virtual, i.e. off mass shell. For 
this purpose the matrix element used in the photoproduction Mj taken at the photon virtuality 
not equal to zero, where i is the index corresponding to the photon polarization, is squared and 
convoluted with the spin averaged electron current squared in the following way: 

I^P = E ^ ^ ^1 M^mM*{Q% (13) 

where \A\'^ is the square of amplitude for the electroproduction of D*-mesons, k\ and /c2 are the 
initial and final momenta of positron, and = — (fci — ^2)^- 

The results of calculations are presented in Fig. 7 for the production of D*-mesons in deep 
inelastic scattering of e~^p in comparison with the experimental data of ZEUS Collaboration 
in the following kinematical region: 1 GeV^ <Q^ < 600 GeV^ -1.5 < r] < 1.5, 1.5 GeV < pr < 
15 GeV, 0.02 <y < 0.7, where y = W^/s. 

Calculating the matrix element we use the running coupling constant of QCD at the scales 



f^R = v"^D* + (the upper curve) and fiR = JAm'jj, + (the lower curve), and the structure 



function of gluon inside the proton is evaluated at the scale \if = \J^i^'d* + Q"^ in the CTEQ4 
parameterization |^. Such the choice allows us, first, to get a good description of electroproduction 
and, second, to reach the scales used in the photoproduction at = 0. 

The ratio of (0(8)) = 1.3 has been supposed the same as in the photoproduction in the 

good agreement with the data. 
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IV. THE TOTAL CROSS SECTION OF CHARM PHOTOPRODUCTION. 



Since we see quite a good description of current data on the photoproduction of charm at 
Pt ~ in the model, it would be of interest to test its predictions for the total cross sections of 
charm photoproduction in a broad region of interaction energy |2^. We evaluate the total cross 
sections of charmed particles at the 7p-energies in the interval from 20 to 300 GeV (see Fig. 8). 

To estimate the total charm yield we multiply the D*^ and D*^ cross section by the factor 
of 1/0.44 being the inverse fraction of these particles in the charm events as was measured in the 
e"''e~ annihilation. Furthermore, to compare with the data at low energies we have to take into 
account the additional contribution of diagrams related to the production of D*-mesons off the 
valence quarks in the subprocess 75^, — (eg) + c D*^ + X. The consideration of such term can 
be performed with no introduction of additional parameters because the normalization of the cross 
section is determined by the same factor of {O'^^^). In contrast to the contribution of 7(7 — > {cq)+c+q 
subprocess, the contribution of 7g„ — * (eg) +c weakly depends on the light quark mass m^. However, 
due to the additional condition (|) relating the values of niq and {O^^^) the implicit dependence 
of charm yield on takes place for the production on the valence quark. So, changing the light 
quark mass from 0.3 to 0.26 GeV the contribution of this subprocess decreases approximately by 
1.5 times. The other terms remain unchanged in this way. 

Thus, varying the mass iriq we can adjust the contribution of the ^q^ — *■ (eg) + e D*^ + X 
subprocess due to the normalization factor of (O^-^-^), and the slow overestimate of experimental 
data can be removed by the small correction of value. 

At the HERA energies the contribution of valence quarks is about several percents, and in the 
framework of present paper we consider a further discussion on this question is inappropriate. 

V. CONCLUSION. 

We have presented the model for the production of charmed mesons in the positron-proton 
scattering. The essential ingredient of consideration is the introduction of fast degrees of freedom 
inside the heavy meson for both the heavy and light quarks, whose differential characteristics are 

15 



calculated in the framework of perturbative QCD at the transverse momenta px > ^c- In the 
model the hadronization is described by means of matrix elements in the form of factorizing the 
hard partonic subprocess and the nonperturbative quark operators corresponding to the transition 
of quark state into the hadron. In this approach, the semilocal duality assumes that the differential 
characteristics of quark pair composing the meson repeat the distributions of hadron in the region 
of applicability of perturbative QCD px > ^c- At large px ^ mc-, as we expect from the common 
theorem on the factorization of hard processes, the cross section has the form of convolution for the 
cross section of heavy quark production with the fragmentation function. 

The region of px ~ is described by additional terms different from the fragmentation. They 
have the asymptotics l/px- We enhance the charm yield in this region due to the contribution 
coming from the color-octet cg-state. 

Within the model uncertainties we have described the data on the photoproduction and elec- 
troproduction of _D*-mesons at the HERA energies as well as the data on the total cross sections of 
charm photoproduction at fixed target. 

The improvement of D*-meson spectrum description in the present model in comparison with 
other approaches is achieved due to the clear physical effect: the interaction of virtual c-quark 
with the fast light quark is increased in the color-octet state by the diagrams of recombination 
in the charm production at px ~ fnc- At large px this contribution becomes negligible and the 
fragmentation regime dominates. It is worth to mention that in our approach the calculation of light 
quark distribution at large px ~ w^c has been done perturbatively. If one neglects the interference 
terms, then this mechanism for the enhancement of charm yield in the forward direction can be 
interpreted as the c-quark interaction with the hadronic remainder (the sea of light quarks) , which 
increase the forward production (so-called " drag-effect" ) |^ . 



The most essential and main advantage of the method used is the possibility to take into account 
the interferences of different contributions as well as the higher twists over the transverse momentum. 

We insist on the use of 'fragmentation' term as not equivalent to 'hadronization' in a general 
case. For example. Fig. 5 shows the most successful theoretical description of rapidity distributions 
for the photoproduction of D* , and the deviations at large rapidities are not so dramatic as in 

16 



other calculations cited above. They explored a point of view coming from a believement: the 
hadronization implies the fragmentation. In order to describe the fragmentation of heavy quark, 
nothing is required beyond the perturbative QCD cross section of single heavy quark convoluted 
with a well-justified function by Peterson et al., say. To consider the general hadronization of fast 
components in the heavy meson, one should introduce some physical quantities defining the fastness 
(cut off low virtualities, or the constituent mass of light quark) and the general probabilities to find 
the fast components inside the meson in the color singlet and octet states. So, the model parameters 
are degenerate into the minimal number (3) of unavoidable quantities in the hadronization of heavy 
meson. Moreover, the data on the electon-positron annihilation into D* fix the sum of octet and 
singlet terms. The normalization of cross section for the photoproduction of D* at large transverse 
momenta and high energies of photon-nucleon collisions restricts the variation of constituent mass, 
while the pT-distribution measures the normalization of octet contribution. Then we reproduce the 
differential cross sections on the rapidity in various ranges of transverse momenta (3 distributions) 
and photon-nucleon collision energies as well as in DIS (6 distributions) with a good accuracy. This 
fact imphes the success of the model in contrast to other QCD-inspired approaches based on the 
fragmentation: they use the variation of charmed quark mass to get a good normalization with 
large diviations in distributions or fit the high px distribution with a loss of correct noramalization 
of integrated cross section. 

It is quite evident that the offered scheme of evaluation for the additional contribution into the 
production of heavy mesons can be transformed to the case of hadronic production. In this sense, 
it would be of interest to compare such estimates with the data on the S-meson yield at the FNAL 
Tevatron, which we plan to make in the forthcoming papers. 
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Fig. 1. The diagram describing the production of D*-meson in the process of 
e+e" J* , Z ^ D* + X to the leading order. 
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Fig. 2. The leading order diagrams for the production of (cg)-state in the 57*- 
inter actions. 
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Fig. 3. The distribution over the transverse momentum of D*-mesons produced 
in the singlet cg-state (the upper sohd curve) in comparison with the 
fragmentation model (the down solid curve) in the 75'-interaction at the 
energy 200 GeV. The analogous distributions are shown for the D*- 
meson production in the octet c^-state (the dashed curves). 
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Fig. 4. The differential distributions in the photoproduction of Z)*-meson over the transverse 

momentum (px) and the pseudorapidity (77) in comparison with the ZEUS data at 

ISO <W < 280 GeV and < 1 GeV^. The dashed hne presents the contribution 

of the color-singlet term, the dotted line does the octet term, and the sum of terms 

is shown by the solid line. The upper and down curves correspond to two scales for 

the factorizing the matrix elements of quark operators as it is described in the text. 
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Fig. 5. The differential distributions in the photoproduction of D*-mesons at 
80 <W < 120 GeV and < 0.01 GeV^ The notations are the same 
as in Fig. 4. 
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Fig. 6. The differential distributions in the photoproduction of Ds and D*- 
mesons at 80 < < 120 GeV and < 0.015 GeV^. The notations 
are the same as in Fig. 4. 
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Fig. 7. The differential cross sections for the production of D*-mesons (a - over 

the photon virtuality in GeV^, b - over the Bjorken x, c - over the 

invariant mass of hadrons in the final state, d - over the transverse 

momentum, e - over the pseudorapidity, / - over the Feynman x) in 

the deep inelastic e'^p-scattering in comparison with the data of ZEUS 

Collaboration. The notations are the same as in Fig. 4. 
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Fig. 8. The total cross section of charm photoproduction versus the energy of 
7p-interaction W. The soUd markers represent the data, and the dashed 
and soUd curves show the model predictions with no valence contribution 
and its addition, respectively. 
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